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The stable periòdic orbits of an area-preserving map on the 2— torus, which is formally a variant of 
the Standard Map, have been shown to explain the quantum accelerator modes that were discovered 
. in experiments with laser-cooled àtoms. We show that their parametric dependence exhibits Arnol'd- 

f**^ ' like tongues and perform a perturbative analysis of such structures. We thus explain the arithmetical 

organisation of the accelerator modes and discuss experimental implications thereof. 
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Experiments of àtoms òptics have discovered the new phenomenon of "quantum accelerator modes" . A subsequently 
formulated theoryQ shows that these modes correspond to stable periòdic orbits of the formally classical dynamical 
[ system, that is defined on the 2-torus by the map: 
OO ■ 

J n +i = J n + 27TÍI + ksin(6 n+ i) mod(27r) , 

e n+x = 6 n + J n mod(27r) . (1) 

This map is a variant of the Standard Map, to which it reduces for íl = 0, and its periòdic orbits will be characterized 
in this paper by two integers p,m so that p is the period and m/p is the winding number "in the J direction". It 
\ should be mentioned that Q does not emerge from the classical limit fi, — > of the atòmic dynamics· and also that the 
quantum accelerator modes are unrelated to the well known accelerator modes of the Standard Map[3,0|, because they 
' do not result of múltiples of 27r being accumulated by an orbit as it winds around the torus in the J direction. In fact 
they also arise of orbits with m = 0, and their origin is subtler; we defer the interested reader to Ref.0- The modes 
reported in correspond to orbits with p = 1; however, the theory predicts that also orbits with higher p should give 
rise to accelerator modes; and such "higher order" modes were indeed observed in subsequent experiments 0. This 
■ opened the way to "accelerator mode spectroscopy" , i. e. systematic classification of modes according to their numbers 
. ', P; m - Then the question arose, which winding ratios m/p correspond to observable modes, and why The answer to 
this question has been recently announcedQ and is presented in full technical detail in the present paper. We show 
that the accelerator modes bcar an analogy to the widely studied mode-locking phenomenon, which is observed in a 
variety of classical mechanical systems 6J. This analogy includes important aspeets such as the Arnol'd tongues, and 
the Farey organisation thereof. At the same time, the present problem has significant differences from well known 
instances of mode-locking in the physical literature, such as, e.g., those which are reducible to the Circle MapQ- 
These differences stem from the fact that Q is a non-dissipative (in fact hamiltonian) dynamical system. 
In this paper these issues are analyzed in detail, thus providing a backbone for the results announced in [2j. We 
develop a perturbation theory for the tongues near their vèrtex, and a heuristic analysis for the "critical region" 
where they break. Based on such results we describe and explain the Farey-like arithmetical regularities that emerge 
from classification of the observed quantum modes, and show that such regularities are encoded by the arithmetical 
process, of constructing suitablc sequences of rational approximants to a real number, which is just the gravity 
acceleration (measured in appropriate units). 

Our perturbative analysis exposes a formal relation to the classical Wannier-Stark problem of a particle subject to a 
constant field plus a sinusoidal field. This relation has quantum mechanical implications, which are discussed in 

This paper consists of two parts. In the lst of these (sections llllllll llVf) we perform analytical and numerical 
analysis of the tongues. Based on these results, in the 2nd part (section[Vj we turn to connections to experiments. 
The most technical aspeets are deferred to Appendices. 
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II. PHASE DIAGRAM. 

The "phase diagram" in Fig^shows the regions of existence of several stable periòdic orbits with different p, m in the 
plane of the parameters fí, k. The origin of stable periòdic orbits of Q associated with any couple p,m of mutually 
prime integers is easily understood. For rational íl = m/p and k — 0, the map has circles of period-p points. As 
generically predicted by the Poincaré-Birkhoff argument Q , at nonzero k these circles are destroyed, and yet an even 
number of period-p points, half of which are stable, survive in their vicinity. At sufhciently small k > such stable 
periòdic orbits exist in whole albeit small intervals of vàlues of fi around m/p. Exactly this fact gives birth to the 
experimcntally observed accelerator modes; indeed, a stable (p,m) orbit of Q with fi in the vicinity of m/p gives 
rise to a quantum accelerator mode, whose physical acceleration is proportional to |íl — m/p| (0; see also sectlVll. 
The persistence of a given winding ratio m/p in a whole region of the space of parameters k, Í7 is where an analogy 
to the "mode locking" may be seen. As shown in Fig^ near the k = axis these regions ("tongues") turn out in 
the shape of wedges, with vertices at fí = m/p, k = 0. The wedges exhibit, at their vèrtex, an angle, and not a cusp, 
as is instead the case, e.g. with the Circle Map, and with systems that reduce to it due to dissipation[l(|. Moving 
to higher k inside a tongue, the periòdic orbit turns unstable, causing the wedge to break and ramify. Bifurcations 
follow, which give rise to swallow-likc structures. Such "critical structures" of different tongues intertwine and overlap 
in complicated ways. A tongue is usually overlapped by others, even before breaking, so stable orbits with different 
p,m coexist; according to numerical computations, such overlaps persist at very small vàlues of fc, marking one more 
diffcrcnce to the usual scenario. It should be noted that higher-period tongues in FigQ]hide lower-period ones, and 
this concurs with graphical and numerical resolution in effacing much of the fine structure of the critical regions. 

III. PERTURBATION THEORY. 

We consider the case when |fc| is small, and is close to a rational number m/p, with m,p mutually prime integers. 
We then write 

ü = m/p + ea (27r) _1 , k = ek . (2) 

where e is a small parameter. The sign of a and e is arbitrary, and k may be assumed nonnegative with no limitation of 
generality. By working out canonical perturbation theory at lst order in e, we determine the finite angle at the vèrtex 
of the p,m tongue, and obtain an estimate for the area of the stable islands. The whole procedure is an adaptation 
of Chirikov's clàssic analysis an d for SI = our results reproduce well-known ones for the Standard Map. 

A. Setup. 

To open the way to a Hamiltonian formulation, we first of all remove mod(27r) from the lst eqn. in and thereby 
translate into a map of the cylinder parametrized by ( J, 9) G M x [0, 2ir) on itself. Doing so, period-p points on 
the torus are turned into non-periodic points on the cylinder, due to the constant drift 27rm/p in the lst eqn. in 
which is not any more suppressed by mod(27r) . For this reason we change variables to L n = J n — 27rrim/p and thus 
obtain: 

L n+ i = L n + ea + ek sin(é> n+ i) , 

9 n+1 = 9 n + L n + 2vrmn/p mod(27r) . (3) 

This defines a map A4 n : (L n , 9 n ) ~~» (L n+ i, # n +i) which explicitly depends on the "time" n. However, A4 n +p = -M. n , 
so, denoting L — L np ,8 — 9 np and L = i( rl +i) p ,é' = é>(, i+1 ) p , the map M 1 -^ : {L,8) (L,9) is defined as in 

M {v) = M np +p-i o M nv+ p-2 o . . . o M np (4) 

and does not depend on n any more. The search for period-p points of Q is thus reduced to search for period-1 
points of A4^ v \ For e — 0, these fill the circles L — R p s , where 

R p , s = n(2 S - X (p))/p , s e Z . (5) 

Here x(-) is the characteristic function of the even integers. We next write (0} at lst order in e in the form of a 
canonical map that affords implementation of canonical perturbation theory. It is easily seen that, at lst order in e, 
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FIG. 1: Phase diagram of the map showing the regions of existence and stability ("tongues") of several periòdic orbits 
with different (p, m). The thick dotted black lines represent the locus of the parameter vàlues used in experiments. The dashed 
black lines show the perturbative theoretical prediction 11411 for the margins of a tongue. 
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the map writes: 



L = L + ea p + efc sm(6 + sL + 7rms(s - l)/p) 



s=l 

= L + eap - efc — G(p,m, 9, L) , 

p — 1 r 

9 = 9 + pL + ea p(p - l)/2 + x(p)tt + efc sin(0 + si + nms(s - l)/p) 

V— 1 S — 1 

= + pL + ea p(p - l)/2 + X (p)tt + efc (— - P^)G(p, m, 9, L) , (6) 



where 



p 

G(p, m, 0, L) = ^{é 6 Q{p, m, L)} , Ç(p, m, i) = ^ e — (-D/p+-í . ( 7 ) 

s=l 

The sums Ç/(p,m, L) are a generalized version of the Gauss sums that are studied in number theory. They play an 
important role in the present problem, and their moduli and arguments will be denoted A(L) and respectively, 
omitting the specification of p and m whenever not strictly necessary. The map (jüj) is not a canonical one, but may 
be turned canonical, at the cost of higher order corrections only, by replacing L by L in the 2nd equation. To show 
this we note that the function 

S(8, L)=0(L~eap)+ X (p)irL+ ipl 2 - ea L ^^~ — + efc G(p, m, L, 6) 

generates a canonical transformation (L, 6) — > (L, 9), given in implicit form by 

- d - 

L = L + ea p - efc — G(p,m, L,6) 
oo 

6 = + pL + X (p)7T - ea p(p + l)/2 + efc ^=G(p, m, L, 9) , (8) 

oL 

provided that the lst equation may be uniquely solved for L. This is indeed the case whenever 

leKlfcrMp^Ml+p/2)]- 1 . (9) 

where c is a numerical constant of order unity. This follows from \t\L — 1| < |efc dQ(p, m, L)/dL\ and from estimate 
(|37fl in Appcndix IVI Cl It is easily seen that replacing L by L in the argument of G in JHJ exactly yields ©■ As G is 
scaled by e, this replacement involves an error of higher order than the lst. 



B. Resonant approximation. 

At lst order in e, the map JHJ may be assumed to describe the evolution associated with the time-dependent, "kicked" 
Hamiltonian 

^ oc 

H{t) = -pL 2 -eap6-aepL/2 + x(p)TrL + ekG(p,m,L,9) ^ S(t - n) . (10) 

n— — oo 

from immediately before one kick to immediately before the next one. This Hamiltonian is a multi-valued function 
on the cylinder, however multi-valuedness disappears on taking derivatives in the Hamilton equations, and so 11UI) 
uniquely determines a "locally Hamiltonian" flow. We change variable to Lq = L — ea /2 and drop inessential 
constants, as well as corrections of higher order in e; and then, in order to remove explícit time dependence, we 
move into an extended phase space with canonical variables (9, tp, Lq, Mq), and therein consider the time-independent 
Floquet Hamiltonian : 

1 x 
H F (0,ip,L o ,M o ) = -pLl + x(p)TrL -eap9 + 2irM Q + ekG{p,m,L a ,9) £ e lm ^ (11) 
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The variable <p is the phase of the periòdic driving, and changes in time according to <p(t) = p(0) + 2nt. In particular, 
eqn.© is obtamed with <^(0) = 0—. We consider Hll|) as a perturbation, scaled by e, of the unperturbed Hamiltonian 

H = l;pLl + x(P)tt£ + 2ttM . 

Points with Lq = Rp. s (ep.©) and arbitrary 9, <p, Mq, are fixed under the evolution generated by Hq in unit time. 
For £^0a 2-parameter family parametrized by ip, Mq survive near Rp, s - These points may be analyzed by Standard 
methods of classical perturbation theory £| in the vicinity of each resonant value i? PiS of the action Lq. This calculation 
is reviewed in Appendix IVI Bl The final result is that, for sufficiently small |e|, and near each resonant action i? PjS , the 
motion in the Lq, 9 space is canonically conjugate at lst order in e to the motion described by the simple Hamiltonian 
in (|12|l below. This result is achieved by three subsequent canonical transformations. The first of these removes the 
oscillating (ip— dependent) part of the perturbation to higher order in e, except for a "resonant" part, by moving 
to appropriate new variables r di,ípx,L\,M\. The 2nd transformation leads to variables ^2,^2,-^2,-^2 such that the 
#2, L2 motion is decoupled from the tp2, M% motion. A final transformation leads to variables #3, L3 such that the lst 
order perturbation term in the Hamiltonian depends on the angle variable #3 alone. The final Hamiltonian is that of 
a pendulum with an added linear potential: 

H rcs = ^pL 2 3 + eV(9 3 ) , V(e 3 ) = -pae 3 + ky/pcos(6 3 ) . (12) 

A previous remark about multi-valuedness of (|1U|) applies to this Hamiltonian, too. In spite of being ill-defined on 
the cylinder, it defines a locally Hamiltonian flow. Replacing the angle 9 by a linear coordinate turns (|12[) into 
the Wannier-Stark (classical) Hamiltonian for a particle moving in a line, under the combined action of a constant 
field and of a sinusoidal static field 0|. Relations between the present problem and the Wannier-Stark problem are 
discussed in @ . 



IV. TONGUES. 



A. Stable fixed points. 

Equilibrium (fixed) points (£3,6*3) of the Hamiltonian l|12fl must satisfy 

L* 3 = , V'iOl) = -pa - fcVpsin(^) = , (13) 

hence they only exist if |a| < fcp -1 / 2 , or, equivalently, 

|í2-m/p| < (2^r l· |fc| p- 1/2 . (14) 

Under strict inequality, i|13|) has two solutions, and one of them is stable. The presence of higher-order corrections, 
(which were dismissed along the way from to Q12[l) turns the dynamics from integrable to quasi-integrable, so, 
assuming a conventional KAM scenario, one may predict a stable orbit of l@J near each resonant action Rp. s , for 
sufficiently small |e|. In order to determino the equilibrium points in the original variables, one may work backwards 
the canonical transformations specified in Appendix IVI Bl and in the end recali L = Lq + ae/2; or else one may directly 
solve for the fixed points of © at the lst order in e. In either case one has to use formulae (|34fl and i|35|) in Appendix 
rvTCl It is thcn found that 

L* = Rp.o + o(e) , (p odd) , 

L* = Rp,o + -ae + -kesin(9*) + o(e) , (p even) , 

9* = - arcsin(aVp/fc) - £(flp,o) + 0(e) . (15) 

The phases ^(i? PjS ) were computed in closed form by number-theoretic means by Hannay and Berrv[Ï4). A chain of 
p fixed points of period p are then obtained for the original map @J on the torus. For small e, these points belong 
to a single primitive periòdic orbit of @J, because they result of a continuous displacement, scaled by e, of points in 
a primitive periòdic orbit of (QJ for e = 0. In the (fi, k) phase diagram, (fT4*)l is satisfied in a region bounded by two 
half-lines originating at k = 0, il — m/p. At small k the half-lincs exccllcntly reproducc the side margins of the (p, m) 
tongue, as determined by numerical calculation of the (p,m) stable periòdic orbits of the exact map @J (see Fig|5], 
and the dashed lines in Fig^J. For p = 1 (|14|) coincides with an exact condition given in 0, which is vàlid at all e. 
For p > 1 it significantly strengthens that condition, but is only vàlid at lst order in e. 
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FIG. 2: Comparison of the map dynamics Q with the resonant Hamiltonian dynamics ll i 2^ in the vicinity of the resonant 
action R p ,o ~ 0, with ea = —0.013, k = 0.1257, p = 5, m = 2. The coordinates shown on the axes are J, 6 for the map dynamics 
(points) and Z/3, #3 — £(0) for the pendulum dynamics (full lines) (cf. eans. l|15|l 'l. 8* and 61 mark the unstable and the stable 
equilibrium points of \\\'1\> respectively; the thick line is the separatrix, and 6* marks its return point. These are marked by full 
circles. 

B. Size of perturbative islands. 

Elliptic motion around a stable equilibrium point generates a stable island in the (£3,03) phase space, and hence 
a stable island for the discrete time motion in the (J, 9) phase space, with the same area (at lst order in e). In 
Fig|21we show a stable island of the exact map (JIJ along with a phase portrait of the Hamiltonian flow (|12fl . In the 
perturbative regime, where the approximation (|12(l is vàlid, we roughly estimate the area A of an island by the area 
enclosed within the separatrix of the integrable pendulum motion (also shown in FigO - T° this end we introduce a 
(positive) parameter A = \a\^/p/k, so that lines A =const. are straight lines through the vèrtex of the (p,m) tongue. 
The axis of the tongue corresponds to A — and the side margins to A = 1, so condition 114(1 is equivalent to < A < 1. 
The estimate is then (see Appendix IVl D|) : 

^«cr 1/4 |fc| 1/2 /(A), (16) 

where c is some adjustable numerical factor of order unity, slowly varying with k and A, and f{\) is an implicit 
function, dcfined in Appcndix IVI Dl the form of which may be inferred from Fig|21 It monotonically decreases from 
8tt to as A increases from to 1, and near these endpoints it behaves like: 

/(A) - 87r-4(47r) 1/2 A 1/2 as A^0 + 

/(A) ~ 3 3 / 2 x 2 7 / 4 (l-A) 5 / 4 as A-»l - . (17) 

Thus, along lines drawn through the vèrtex of a tongue, and sufficiently close to the vèrtex, A decreases proportional 

to \fk. The estimate l|16[) and the asymptotics l|17[) are well confirmed by direct numerical estimation of areas of 
stable islands of , as shown in Figs [21 El E] 
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C. Limits of validity. 

A crude upper bound for the validity of perturbative analysis is set by overlapping between islands, belonging in 
the same mode and in neighbouring modes as well. If only the former type of overlapping is considered, then the 
no-overlap condition reads pA < 4n 2 and yields |fc| < const.p~ 3 / 2 . Turning estimates based on the overlapping 
criterion into exact (albeit possibly non-optimal) ones is quite problemàtic 4]. However, one may assume that the 
dependence on the period p is essentially correct. Two further conditions are set by the validity of © itself as a lst 
order approximation to which results in the bound Q, and by the validity of the resonant approximation, which 
results in the bound (|31|) . The logarithmic corrections in i|üll. (|o4[) are likely to be artifacts of our derivation; in any 
case, both bounds have nearly the same dependence on the period p as predicted by the "overlapping criterion" . 

D. Crisis of the tongues. 

The perturbative estimate l|16(l is vàlid near the vèrtex of a tongue. On further moving upwards in the phase 
diagram, the area of stable islands first increases up to a maximal value, and then it decreases through strong 
oscillations (FigQJ. The islands finally disappear, as soon as an upper critica! border of stability is reached (Fig^ 
Fig|SJ). On trespassing this border bifurcations are observed 0, giving rise to stable, primitive (p, m) orbits, with p, m 
non relatively prime. The morphology of tongues in the critical regions where they break is superficially remindful 
of that observed with other maps[Ï5j and its analysis is outside the scope of this work. In the case p = 1, exact 
non-perturbative calculation of the fixed points and of their stability is possible, showing that the upper stability 
condition involves the 2nd order in e 0] . Stability thresholds estimated from the trace of the derivative of the map 
© at the fixed points miss effects of higher order corrections that were neglected in deriving (jSJ) itself from QJ. 
We therefore resort to numèrics. Having in mind the border (|14|l and the discussion in sec II V Cl we refer each 
tongue to scaled variables |fc|p 3 / 2 , \fl — m/p|p 2 . The horizontal scaling is chosen such that all tongues have the 
same vèrtex, and the same angle at their vèrtex. Figjü] shows that the subcritical parts of all inspected tongues 
occupy roughly the same region in the plane of the scaled variables. A similar indication is given by Fig0| It is 
worth noting that scaling with the variable p 3//2 |fc| is predicted by l|16|) for the total area pA of the islands of a 
period-p orbit, in the perturbative regime of small k. On the basis of all such indications we assume that the critical 
region where a tongue breaks is roughly located around |fc C R.| — b p~ 3 / 2 , with ò ~ 6 as suggested by Figs l5l4l The 
critical border that way (somewhat vaguely) defined has the same dependence on p as the previously discussed borders. 




FIG. 3: Numerically estimated area of the stable island of the (1, 1) orbit of the map at constant k = 0.3666, as a function 
of the distance from the centre of the tongue, as measured by the parameter X = |a|p 1//2 /fc . The full line shows the theoretical 
estimate Il·l^ , with c = 0.66. Inset: bilogarithmic magnification of the shaded region near X = 1. The slope of the dashed line 
is 5/4, as predicted by 1171 . 
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FIG. 4: Behaviour of the numerically estimated total area of the chain of stable islands of a periòdic orbit, as k increases 
along the line A = 0.1405 (dotted line in Fig.5). Results are shown for several orbits with periods p < 18. Inset: bilogarithmic 
magnifïcation of the small-fc shaded region. The full line is the perturbative prediction from (llbH with c = 0.68. 



V. SPECTROSCOPY OF TONGUES. 

The theory developed in the previous sections provides a quantitative description of the gross structure of the phase 
diagram. Application to quantum accelerator modes is made in this Section. Wc first elucidatc the physical meaning 
of the phase diagram in this context. 

A. Experiments with cold àtoms. 

In experiments pl |2Ï| . cold caesium àtoms of mass m are subject to very short pulses, or "kicks", with a period T 
in time. The strength of a kick periodically depends on the position of an atom (assumed to move in a line) at the 
kicking time, with a period 2tt/G . Its maximal value is denoted k. Inbetween kicks, an atom freely falls with the 
gravitational acceleration g. The accelerator modes are observed when T is close to special resonant vàlues, which 
are given by Ti = 2ním/ {HG 2 ) with i any integer. Writing T = 2^(1 + e/(27r£)), the small parameter |e| is found 
to play the formal role of a Planck constant in the quantum equations of motion 2] . In the limit when this Planck 
constant tends to the atòmic dynamics are governed by the "e-classical" [2^ map with 2nQ — GT 2 g, k = fce , 
and J = ne, where n is the atòmic momentum measured in units of hG |23|. The theory shows that àtoms which are 
trapped in a stable island of the map move with constant physical acceleration, thereby giving rise to an accelerator 
mode. Their acceleration relative to that of freely falling àtoms is given, in units of h 2 G 3 /m 2 , by the parameter a 
|24| . The acceleration a of a mode may be inferred from the experimental momentum distributions of the àtoms after 
a given number of pulses. As is known, the rational winding number r — m/p is then determined. The integers p 
and j =sgn(e)m have been respectively termed the order and the jumping index |2j| of a mode 0- 

B. Quantum Phase Diagram. 

At nonzero vàlues of the "Planck constant" e, the e-classical picture is subject to quantal modifications. Whilc the 
e-classical dynamics only depend on two parameters k, íl, the quantum dynamics additionally depend on the "Planck's 
constant" e, which is not determined by k and Q alonc. Thus, for instance, the acceleration a of a mode is not a 
e-classical variable, because its value at any given point (k, Q) in a tongue further depends on e, which is a priori 
arbitrary. Once a value is chosen for e, the horizontal width of the (p,m) tongue at any given value of k, multiplicd 
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FIG. 5: The lower, subcritical parts of tongues of different periods p < 29, once rescaled as indicated on the axes, reduce into 
roughly the same region. The dashed lines show the analytical stability borders of the (1, 0) tongue. The dotted line is drawn 
for reference to Fig0] 



on 7r/|e|, yiclds the maximal (in absolutc valuc) acceleration that may be attained in the (p,tn) accelerator mode with 
the given k. 

Quantum effects would efface fine structures in the phasc diagram, if determined by too small islands compared to 
|e|. Hence, if a value of e ^ were chosen once and for all, independently of the vàlues of k,íl, then high-period 
tongues would be quantally irrelevant, and low-order modes might be observed only in the inner parts of tongues, 
sumcicntly far from their borders, where the islands shrink to zero. However, in experiments, e is not fixed, as k is 
varied by changing e at constant k. As shown by the estimate (|16fl . in this way the area A of an island decreases 
with yfef, so the ratio A/\e\ grows arbitrarily large at small k. Consequently, the e-classical dynamics are more and 
more accurately reflected in the quantum dynamics of àtoms, as the vèrtex of a tongue is approached. In particular, 
quantum effects do not set restrictions of principle to the observation of modes of arbitrarily large order. On the 
contrary, the breakdown of a tongue occurs relatively far from its vèrtex, and quantum effects may not be negligible 
there. Significant quantal modifications on the e-classical critical behaviour have been observed and discussed in 0- 



C. Arithmetics of accelerator modes. 

A genèric feature of mode-locking phenomena is classification of the locked modes by means of the Farey tree0]> which 
is a Standard technique in number theory^ïj- This construction is based on a curious property^^ of the irreducible 
fractions and orders the rational numbers in a hierarchy, which turns out to essentially reproduce the natural hierarchy 
of the modes, as dictated by their physical "visibility" . In this subsection we discuss the arithmctical organisation of 
the quantum accelerator modes. To fix ideas, we assume that Í2 ranges in the interval [0, 1]. For r a rational number 
in [0, 1] we denote m(r)/p(r) the corresponding irreducible fraction. 
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1. Farey rule. 

Due to finiteness of the interaction time, only modes up to some maximal order M can be detected in experiments 
and in numerical simulations; the integer M being roughly determined by the interaction time |26j . The set Tm of 
all rational numbers r such that p(r) < M, arranged in increasing order, is known as the M-th Farey series. If a 
rational r is thought of as the winding ratio of an orbit, then Tm provides a catalog of all the tongues of period < M, 
ordered from left to right in the phase diagram. This may be termed the M-th Farey family of tongues. If M > 1, 
then statements (FI) and (F2) below are true of the Farey series Tm US : 

(FI) If r*i and r-i with r-i > r% are nearest neighbours in the Farey series Tm, then p(ri)m(r2) — p(?~2)tTi(?"i) = 1. In 
particular, p(ri) and p^) are relativcly prime integers, and so are m(ri) and m(r2). 

(F2) If f3 is the element of Tm that follows r-i on the right, then r^ is the Farey mediant of r% and ra, i.e.: 

_ m(ri) +m(r 3 ) 
P(ri)+P(r 3 ) 

These facts have the following consequence. Let two tongues r\,r% be "adjacent", in the sense that no other tongue 
exist between them, with a period less or equal to the largest of the periods p(ri),p(r2). Then (F2) implies that the 
tongue with the smallcst period, to be found between them, is the tongue associatcd with the Farey mediant of r\ and 
r2 . This may be called the "Farey rule" . A qualitative formulation of this rule is : the next most visible tongue to be 
observed between two adjacent tongues labeled by rationals r\ and r2 is the tongue labeled by the Farey mediant of 
n and T%. This rule does no more than reflect the fact that tongues are the less visible, the higher their period. 



2. Experimental Paths. 

Orders and jumping indexes of quantum accelerator modes are identified, by monitoring the atòmic momentum 
distributions that are obtained after a fixed time with different parameter vàlues. The latter are generated by 
continuously varying a single control parameter: the pulse period T (or, equivalcntly, e) in ranges close to the 
resonant vàlues Ti. The locus of the corresponding points in the (fc, Í2) plane is then a continuous curve, that will be 
termed the Experimental Path (EP) in the following. The problem then arises, of classifying the tongues, that have 
a significant intersection with a given EP. These are but a small subclass of the Farey family of tongues that fits the 
given interaction time, because many tongues in that family are met by the EP in their overcritical regions, where 
islands are typically small, yielding hardly if at all detectable modes. Thus the analysis of observable modes rests on 
three key facts. The first of these is existence of a critical border fsect lIVD|) . Second, an EP hits the k = axis at a 
value of Q, given by u> = GgTf / {2it), which corresponds to e = 0, or T = T(. Third, the quantum dynamics grow more 
and more quasi-e-classical as k — is approached along an EP, and this justifics analysis based on the e-classical phase 
diagram. These three facts reduce the problem of predicting the observable modes to a number-theoretic problem, of 
constructing suitable sequences of rational approximants to the real number u>. If space is measured in units of the 
spatial period of the kicks, and time in units of Tt, then oj is just the gravity acceleration. 

An EP may in general be described by an equation Çl = lj + <í>(fc), where $(fc) is some strictly monotonic function 
such that $(0) = 0. The experiments in |5j will be used as a model case in this section, and the corresponding EPs 
are shown by the black dotted lines in Fig^. Each choice of l = 1, 2, 3 yields an EP consisting of two lines which will 
be henceforth referred to as the two "arms" of the EP. The left (resp., right) arms of the EPs corresponds to negative 
(resp., positive) vàlues of e. With 1 = 2 the arms meet at k — 0, SI = u> w 0.390152... and are approximately linear at 
small k: \k\ ps a\u> — fi|, with a — h 2 G 3 k/(2m 2 £g); so we may assume $(fc) « a k at small vàlues of e for the case 
of Figffl 

Independently of the specific form of í", the intersection of an EP with the subcritical (p, m) tongue is defined by two 
conditions, one dictated by i|14|) . and the other by the critical border \k CIit \ ~ òp~ 3 / 2 with b ~ 2ir: 

|m/p-íï| < (2^)- 1 p- 1 / 2 |fc| , |w-íï| = |*(À)| < |$(27rp- 3 / 2 )| (18) 

These inequalities lead to 

\u - m/p| < {2n)- 1 bp- 2 + |$(2^p- 3 / 2 )| ~ p" 2 + |$(27rp- 3 / 2 )| , (19) 

and show that the winding ratios r = m/p that are observed along an EP have to approximate lo the better, the 
smaller e. 
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3. Farey algorithm. 

According to (|19f) . the winding ratios of modes oberved along an EP form a sequence of rational approximants to 
lo. This fact was already noted in 2], and the question arose, which of the densely many rational winding ratios that 
lie arbitrarily close to lo are actually observable. The Farey rule may be used to answer this question j^. Modes are 
the more visible, the smaller their order; therefore, modes observed on moving along an EP towards lo should achieve 
better and better approximations to w, at the least possible cost in terms of their orders p. Issues of criticality, and 
of quasi-e-classicality, additionally suggest that, as a thumb rule, modes should be more safely observable near the 
vèrtex of their tongues. These indications suggest the following construction. Assuming that fi ranges in the interval 
[0, 1], the strongest modes are the period-1 ones, (1, 0) and (1, 1). According to the Farey rule in the end of sect IV Cll 
the next strongest mode is associated with their Farey mediant, (2, 1) = (1 + 1,0+ 1). At the next step two furthcr 
Farey mediants (3, 1) = (1 + 2, + 1) and (3, 2) = (1 + 2, 1 + 1) appear. The former is closer to lo ss 0.39... than the 
latter, so its tongue intersects the EP at lower |e| . It is therefore expected to produce a stronger mode; so we discard 
(3, 2), and restrict to the interval [1/3, 1/2]. The process may then be iterated. At the n-th step, it will have singled 
out two rationals, r' n u and w , such that the Farey interval F u n = [r' n ul , r" u ] contains lo, but does not contain any 
rational with a divisor smaller than p{r' nuJ ) + p(r' 7 ' luJ ). If lo is itsclf a rational, then it is eventually obtained as the 
Farey mediant of the endpoints of a Farey interval íLj.n, and then the process terminates. The process just described 
is an arithmetic recursion for generating rational approximants to lo, that will be termed here "the Farey algorithm". 
Out of all possible rational approximants to w, it selects the endpoints of the Farey intervals of lo, as the winding 
ratios whose prediction is safer. 

4- Accelerator modes, as Rational Approximants of Gravity. 

It is now necessary to discuss consistency of the Farey algorithm with the key condition l|19f) . that dictates the rate 
at which modes of increasing order p have to approximate the gravity acceleration lo. This rate depends on the form 
of the function <ï>; however, in no case it is required to be faster than quadratic, thanks to the lst term on the rhs of 
JTÏÏJ). For this reason, observation of the principal convergents (or simply the convergents) to lo is always expected. 
These are the rationals s WjTl that are obtained by truncating the continued fraction of lo, and are "best rational 
approximants" to lo in the sense that (Thm.182 in [Ï9|. p.151): 

p(so,,„)|w - s U!n \ < p(r)|w - r\ , whenever p(r) < p(s ai>ri ) . (20) 

They are known to satisfy \lo — s w ,„| < p(s Wi „)~ 2 , and hence (|19f) as well, and are in fact clearly detected in experiments 
and in numerical simulations. The Farey algorithm generates all of the convergents to lo. As shown in Appendix 
IVI Al at least one endpoint of each Farey interval is a convergent to lo; however, except for quite particular choices 
of lo, the Farey algorithm generates more approximants than just the convergents, and so a Farey interval generated 
by the algorithm may have only one endpoint given by a convergent. In that case, that very convergent is also an 
endpoint of the next generated interval, and possibly of subsequently generated ones, until the construction produces 
the next convergent at the other endpoint. By construction, r' n u is the rational that yields the best approximation 
from the left in the class of all rationals r with p(r) < p(r' n w ); and r' r [ has the the same property in what concerns 
approximations from the right. The one of the two approximants r^ ^ , r" ^ that lies closer to lo is called the n-th 
Farey approximant psl ] to lo, and will be denoted r* n . It is by construction a best rational approximant to lo, in a 
weaker sense than <|20ll : notably, 

|w - <J < k - r| whenever p(r) < p(r* >n ) . (21) 

The approximation of lo which is granted by a Farey approximant may not be quadratic when the approximant is not 
a convergent; therefore, whether a Farey approximant satisfies condition l|19fl depends on the form of the function $. 
That condition is the more exacting, the steeper the EP, which is the graph of the function <E>. For instance, in the 
extreme case when the EP is a vertical line drawn through lo (which corresponds to a — oo in our model case p9| 
), the 2nd term in l|19|l is absent and the approximation has to be strictly quadratic. This may rule out some of the 
non-principal approximants produced by the algorithm, depending on arithmetical propcrtics of lo. For lo equal to 
the Golden Ratio, all the Farey endpoints are principal convergents. All the corresponding modes, and none other, 
were observed in numerical simulations of the atòmic dynamics (Fig|üjl . For a "less irrational" choice lo = ir — 3, 
the Farey algorithm generates many other approximants besides the convergents. All those with order p < 106 were 
detected by our numerical simulations of the atòmic dynamics up to 800 pulses, except for two, and these were found 
to violate l(TÜ)) . 

When the EP is not vertical the 2nd term in (|19fl opens the way to non-quadratic approximants. For the EP we 
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consider here, this term is given by 2na~ 1 p~ 3 ^ 2 and prevails over the lst term at sufficiently large p. For this reason, 
according to theorem (T6) in Appendix IVl Al for almost all u> (in the sense of Lebesgue measure), all of the Farey 
endpoints, except possibly for finitely many exceptions, satisfy (|19|) . Since the 2nd term in (|19|l prevails over the 
lst the later (that is, at larger p), the steeper the EP, (in our model case, the larger a), the possible "finitely many 
exceptions" may be relevant in the analysis of data, which of necessity cannot extend to arbitrarily large p. 
In conclusion, the Farey-based prediction of modes may suffer exceptions, both in the case when the EP is very steep, 
and in the opposite case when it is very flat, by generating more approximants (in the former case), and less (in the 
latter case) than allowed by l|19fl . In addition, l|19|) itself is not exact, as it rests on a coarsely defined critical border. 
In particular, large fragments of broken tongues, too, may produce significant modes at times. 

5. Final remarks. 

Left- and Right-hand Modes. By construction, r u n and r'ú n approximate lo from the left and from the right respec- 
tively. Therefore, the vèrtex of the n tongue lies on the left of ui, so the intersection of the tongue with the left arm 

of the EP lies at significantly lower k than its intersection with the right arm; hence, it should be preferably observed 
at e < 0. It is in fact an empirical observation, that left (resp., right) approximants preferably occur at negative 
(resp., positive) vàlues of e. In particular, two successive convergents to u> approximate uj from opposite sides, so the 
corresponding modes are in principle expected on opposite arms. This is not a strict rule, as tongues with relatively 
small p may have significant intersection with both arms of an EP. Some low-period modes could indeed be observed 
on both arms, and were found to correspond to convergents of lo. However, this cannot happen when the period of a 
tongue is so large that the slope of its margins is larger than the slope of the arm lying on the opposite side 

with respect to ui. Thus, modes of sufficiently large order should never be expected on both sides. In that case the 
two arms coincide, so there is complete symmetry between e > and e < 0, as can be seen in Fig. [ü] 
The case of rational lo. If u> is a rational number r/s , then the Farey algorithm eventually terminates. It is therefore 
expected that, whatever the observation time, only a finite number of modes are observable. This case has been 
experimentally realized, too [2íÍj . The arms of the EP exactly meet at the vèrtex of the (s, r) tongue. If their slope is 
larger than 2tï^/s (as in |20|) then the (s,r) mode is observed on both arms, and at all vàlues of |e| below a certain 
value determined by the critical border of the (s, r) tongue. On the contrary, the (s, r) mode could never be observed 
if 2-Ky/s were larger than the slope of the arms. 
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VI. APPENDIX. 
A. The Farey Algorithm. 

What in this paper is termed "the Farey algorithm" is a recursive means of constructing rational approximants to a 
given real number u>, by iterated calculation of Farey mediants. Though the role of the Farey properties (F1),(F2) 
ísect IV (JÏJ) in the process of rational approximation is a bàsic notion in the theory of numbersjïj]], it is not easy to 
locate references wherein those aspeets which are directly used in this paper be presented in a self-contained way. 
Such a self-contained presentation is given in this Appendix. The one input we use is the Basic Theorem stated 
below, which is equivalent to properties (FI) and (F2) of the Farey series in scct IV CU Proofs of (FI) and (F2), 
hence of the bàsic theorem, may be found, e.g., in ref.[ïj|. 

Let uj 6 [0,1] be fixed. Given a rational r £ [0,1] we denote m(r)/p(r) the corresponding irreducible fraction. We 
further denote d u {r) = \r — lu\ and <5 w (r) = p(r)d w (r). We say that a rational r is S-closer to uj than another rational 
s if Suj(r) < <5 w (s); and we say that r is d-closer to lo than s if d w (r) < d (jj (s). 

DO:j4 best rational d-approximant (dBA) to lo is a rational r such that d u (r) < d ÜJ (s) for any rational s with 
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FIG. 6: A Fibonacci sequence of quantum accelerator modes, as revealed by numerically computed momentum distributions of 
a cloud of 50 àtoms after 400 pulses, with k — 0.8-7T, Q, =const.= u — (VE — l)/2 (the Golden Mean), and for different vàlues of 
the kicking period T = Ti (1 + e/ (27rí)) in the vicinity of T = Tg.,1 — 2. The atòmic momentum n is shown on the vertical axis. 
It is measured in a free-fallin g frame, and is assumed positive in the direction of gravity. On the horizontal axis e is replaced 
by the variable e* = sgn(e)i^e|/27r, which affords relative magnification of the small-e region. The hyperbola-like zones of 
enhanced population are the accelerator modes. Their shapes and positions yield access to their winding ratios m/p : (A) 2/3, 
(B) 3/5, (C) 5/8, (D) 8/13, (E) 13/21, (F) 21/34, (G) 34/55, (H) 55/89, (/) 89/144. These are the principal convergents to 
the Golden Mean. Note the exact symmetry under reflection in the origin of the axes. 

p(s) < p(r). Replacing d by S yields the definition of a best rational 5-approximant (8BA) to lo. 

The dBAs and the óBAs to u> are respectively known as the Farey approximants and the principal convergents to fi. 
From the definition it follows that every SBA is at once a dBA. Another immediate consequence is: 
TO: Let ri and r-i be two successive dBAs to lo, in the sense that p(r2) > p(»r), an d n o rational r with 
p(ri) < p(r) < p(r2) is a dBA. Then d u (r) > d u (ri) for all rational r with p(ri) < p(r) < p(r2)- The statement 
remains true on replacing d by S. 

Dl: A Farey interval is a closed interval with rational endpoints r',r" (r 1 < r" ) satisfying m(r")p(r') — m(r')p(r") — 
1, or, equivalently, r" — r' = l/(p(r')p(r")). 

D2: The Farey mediant of two rationals r', r" is the rational r' r" = (m(r') + m(r"))/(p(r') + p(r")). 

BT (The Basic Theorem): The following statements are equivalent: (a) F — [r J , r"] is a Farey interval, (b) p(r) > 
p(r') + p(r") holds true of any rational r with r' < r < r" ; equality holding if, and only if r — r' © r" . 

Denote the family of all the intervals F Ç [0, 1] such that F is a Farey interval and u> is an interna! point of F. 
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Further define a family of intervals T u so that = T w for irrational ui and T^, — U {[ui]} for rational ui, where 
[w] is a shorthand notation for the interval [w,o>]. 

If F = [r' , r"\ is a Farey interval, then each of the intervals [r' , r' © r"], [r' © r", r"] is a Farey interval. One may then 
define a map 3^ : JF W — * as follows. If F — [r 1 , r"] <E and r 1 © r" ^ uj, then ^(F) is the one of the intervals 
[r'.r'ffir"], [r'©r",r"] that contains w. If F= [r',r"\ e and r'ffir" = w, then&^F) = [co]. Finally fo, ([w]) =Jw]. 
The following proposition shows that the map 3w provides an algorithm for recursively generating the whole of T u . 

Tl (The Farey algorithm) :Let u> G (0, 1) be given. For integer n > define [r' un ,r^ n ] = F^, n = #"([0, 1]). Then: 

(a) {-FL in }n>o is a monotone nonincreasing sequence (in the set theoretical sense) of closed intervals; moreover 
lini,woo \F u , n \ = and n n >o-FL,n = M, 

(b) If lo is rational then F UjTl — [uj] eventually, 

( c) {-Fw,n}n>0 = ■ 

Proof: (a) immediately follows from Dl, D2, and from the definition of 3"^. (b): if F^.n M then F w , n 1S a Farey 
interval and contains ui as an infernal point. Duc to BT(b), the family of such Farey intervals is finite whenever u> is 
rational. 

(c): we have to show that F £ implies F — F u<n for some integer n. If F = [0, 1] then F = F U: q, and if F = [uj] 
then uj is rational and the claim follows from (b). Thus we assume F = [r',r"] C [0, 1] with r' < r", and then (a) 
implies that F C F u n can hold only for finitely many vàlues of n. Let N be the largest such value. From F ^ [uj] 
it follows that F Wj tv ^ M, so r'^ N r" N is an endpoint of F Uí n+i, and then r' u N © r" N g F because F is not 
strictly a subset of F w ^+i by the definition of N . If r'^ N (Br"^ N were an infernal point of F, then BT(b) would imply 
P( r L n © r 'L n) — P( r L n) + P( r " n) — P( r ') + P( r ") which is impossible because one at least of r' and r" is an infernal 
point of F^^n and so its divisor is not less than p(r' u N © r'ú N ). Therefore, r' u N © r'ú N is an endpoint of both F and 
-FL.ív+i· Since the former is not strictly a subset of the latter, and both contain uj, F — i^jv+i follows. □ 

T2:At least one endpoint of each F^^ is a dBA to uj; and every dBA to uj is an endpoint of some F^^. 

Proof: Denote r* the endpoint of F u n that is cí-closer to ui. No rational with a divisor less than p(r*) lies inside F un 

by construction, so r* is a dBA. 

Conversely, let r be a dBA. The claim is obviously true if r = or r = 1, so let r lie strictly inside Fq = [0,1], 
and let m be the largest integer such that r is an infernal point of F Uim . Due to BT(b), m is a finite number, and 
P(r) > PK, m ) + PK,J = P( r L, m © <U)- If P( r ) > P( r 'u.,m © <U)> then <f w (r) < cL(?Í, m © CJ, because r is 
a dBA; hence r is an infernal point of ÍL )m+ i, contrary to the definition of m. Therefore, p(r) = p(r^ ro © r" m ), 
leading to r = m © r" m . Hence r is an endpoint of F m+ i. □ 

(T2) in particular implies, that all principal convergents to uj are generated by the Farey algorithm. The way this is 
done is clarified by (T5) below. The following propositions (T3) and (T4) are lemmata to proposition (T5). 
T3:7f àuir'w n ) = ^oí( r ",n) then uj is rational and F u , n +i = [uj]. If Süj(r' u n ) ^ Sw(r'ú n ) then the endpoint of ' F u>n that 
is 5-closer to uj is also an endpoint of F Utn +%. 

Proof: Fu^n is either [uj] or a Farey interval. In the former case ^(r^ n ) — óuir'ú n ) = and the claim is obvious. In 
the latter case ^ , ò^r'^ n ) ^ 0, and 



which may be rewritten as: 



If &u>(r'u,, n ) = <5 w (r^„), then 



PK,n)PK,n) ' 

<U< JP«n) + *»K, n W U ,n) = 1 ■ ( 22 ) 



w r uj,n — d u [r u „) — r n © r w 

P('t,n)(P(^,n) +P(^,«)) 



7* 



hence r' u n © r" ra = o; and by definition = [uj]. Threfore, if F w , n+ i ^ [ui] then one of the endpoints of F Wj „ is 

ó-closer to ui than the other endpoint. Denoting r* this endpoint, 122(1 implies 

1 

> min {Suir'J^uir")} = p{r*)d w (r^) , 



PK,n)+PK,n) 

We are thus led to 



P(^)(Pfò,J+PW,n)) 
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As r L n © r 'úi n IS an endpoint of F W: „ + i but not of Fu, n the claim is proven. □ 

T4: Let one, but not both, of the endpoints of F Uin be a principal convergent to ui. Then the same principal convergent 
is also an endpoint of F uj n+1 whenever F u n +i ^ [lo]. 

Proof: without loss of generality assume that r'^ n is a principal convergent and that r" n is not a principal convergent. 
The assumptions enforce F u>n ^ [lo]. If S u (r'^ n ) = ^(r^„) then F^ : „ + i = [lo] duc to (T3). If S u (r'ú >n ) > <LK,, n ) 
then the claim follows from (T3). Let us show that ó~ u (r'ú „) < n ) is impossible. If <L(V" n ) < ^('"Lra)) 

then, due to (TO), there must be a principal convergent s with p(r^ n ) < p(s) < p(r" n ), and, due to (T2), s is an 
endpoint of some Farey interval Now F u>m Ç is excluded, because r" n is not a principal convergent, 

and p(s) < p(r'ú <n )- Therefore, F W; „ c F Utm ; but then p{r' un ) < p(s) and BT(b) enforce F u , m = [r^ n ,s], whence 
r L.m+i > r L n because of (T3). Together with F w , n C F u , m +i, this leads to the contradiction r'^ n > m+1 > n . □ 



T5: 4í feasí one endpoint of each PL,n *s a principal convergent to lo . 

Proof: the claim is true of Fu,o- Assume it is true of F w<n . Without loss of generality suppose that r' w n is a principal 
convergent. One of the following is true: 

- -F^n+i = [lo]. Then lo is rational, hence a principal convergent to itself. 

- i^, n +i 7^ [ui], and r" n is a principal convergent, too. The claim follows because F u n +\ has an endpoint in common 
with Fu n by construction. 

- F Uin +i ^ [lo], and r" n is not a principal convergent. Then n — r'^ n+1 due to (T4).D 



It is well known that the principal convergents to an irrational lo provide a "quadratic" approximation to lo, and that, 
for almost all irrationals (in the sense of Lebesgue measure), faster-than-quadratic approximation is impossible. Our 
final proposition states that, for almost all irrational lo, all of the approximants generated by the Farey algorithm 
provide a "quasi-quadratic" approximation at worst. 
T6: For any < n < 1, 

lim p{r'^ n f^d u {r' u J = and lim p{r'^ n ) 2 ' v d^{r^ n ) = for (Lebesgue) almost all lo e (0, 1) . 

n — *oc ' ' n — >oo ' ' 

Proof: consider the lst equality; the argument for the 2nd is identical. Let C C [0, 1] be the set of all the lo such 
that the equality is not true. If lo G C, then 

L(ui) ETkp^J 2 -'^^) > ; 

n — >oo 

so, denoting L'(u>) = 1 if L(u>) — oo and L'(u>) = L(u>)/2 otherwise, the inequality d u (r' u n ) > L' (ui)/p(r' UJ „) 2 ~'' holds 
true for infinitely many vàlues of n. For all such n's, 

l =\ F I > d (r > ) > 

entails p{r' u n ) > (L' (u>)p(r" n ))~ , and hence 

1 1 

) < \F u , n \ = - -. < 



II \2+r,' 



PK,n)PK,n) L"(L0) P (rZ, n ) 

where L"(lo) = (L'(ui))~ and r/ = n/(l — rj) > 0. Hence C C Un>iBn whcre Bn is the set of all lo g [0, 1] such 
that the inequality \ui — r\ < Np(r)~ 2 ~ v with n > holds truc for infinitely many rationals r. Each Bn is known to 
have zero Lebesgue measure. □ 



B. Derivation of the resonant Hamiltonian. 



Let a canonical transformation be generated by a function S — BL\ + ipM\ + eSi(6,Lp,Lx,M\). In order to totally 
remove the oscillating part of the Hamiltonian JHJ, S% ought to solve the equation: 

u (L 1 )^+2n^=-kG(p,m,L 1 ,e) £ e im * . (23) 

' m — — oo 

where uj{L\) = dHo/dL\z J =L 1 — pL± + xip) 71 - Writing the solution as 

S 1 (0,<p,Li,M 1 ) = k Y, ^r, m {Lx,M 1 )e i( - re+m ^ (24) 

r,mÇZ 
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where 



leads to 

*T,m(Li,Mi) = - — l ± , 25 

27TTO + r(pLi + x(P)tt) 

G r (Li) = 7^/ dee- ire G(p,m,L 1 ,6) = l{6 r!l g(p,m,L 1 ) + S r ^ 1 g*(p,m,L 1 )} , (26) 
2tt Jo 2 

Ean. (|25|l cannot be solved if L\ — Rp, s , the resonant vàlues defined in eq.JSJ), and (m, r) = (±s, =Fl). Therefore, in the 
vicinity of i?p, s , only terms with m ^ ±s, r ^ =pl will be removed to higher order, by means of the generating function 
that is obtained by summing (|24J) over (r, m) ^ (=pl,±s) with ay,m given by (|25|1 . Using the Fourier expansioni 

7r * — ' a + m 

which is vàlid for any real noninteger a, this calculation yields the function : 

^ 1 ; 2sin(p(L!- flp,,)/2) 

, . , , sin(6* — + £(íi)) , „ 

+ ' kMM • (27) 

As expected, the transformation generated by this function is singular at L\ = Rp, S ' for s' ^ s; furthermore, it is 
discontinuous at ip = 0, due to the singular nature of the periòdic driving. The "resonant" terms (m, r) = (±s,=Fl) 
remain at the lst order, and sum up to 

6 -^{G(p, m, L x y^- S ^ + G* (p, m, LOe-^-^)} = ek G(p, m, Í -s^L,). 
In this way the resonant Hamiltonian is found, that describes the motion near R p jS at lst order in e: 

fíir res , s = -pL\ - ea p0 x + nx(p)Li + 2ttM 1 + ek G(p, m, 6 X - 8<pi,L x ) . 
One further canonical change of variables : 

di -> 02 = 01- síp 2 , 

(fil — > <^2 = ) 

Mi -> M 2 = Mi + sL 2 

Li — » Í2 = í'l — ííp.s (28) 
decouples the (Í2,fe) motion from the (M2,y>2) motion, and the ^2,62) Hamiltonian reads: 

H res , s = ^pLj - epa8 2 + ek G{p, m, 2 , L 2 + Rp, s ) . (29) 

The Í2— dependence in the 3rd term may be removed to 2nd order in e by one final canonical transformation to 
variables $3, Í3. This is defined by the generating function: 

S 3 (&2, L3) = (9 2 + Ç(R P , S ))L 3 - ek p- 1 ^ A S (L 3 )} , 

where 

A S (L 3 ) = L- 1 [g(p, m, L 3 + Rp,.) - G(p, m, R PtS )} . 

Then, formally, 

L 2 = Ls-efcp-^le^A^La)} , 

63 = 2 + Ç(R P , S ) - ek p^S {e ie *dA s (L 3 )/dL 3 } , (30) 
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Replacing in l|29(l . and dropping inessential constants, one obtains 

1 



, , ... - -pL 2 -ea6 3 + ek A{R PiS ) cos(<9 3 ) + 0{e 2 ) 



which, using (|3(j|l in Appendix IVI Cl yields the hamiltonian in eqn. l(T2|) in the main text. The formal transformation 
(|30|1 is justified provided e is suíhciently small, notably 

líK^ifcp^p^ina+p^)]- 1 , (31) 

where C4 is a numerical constant of order unity. In fact , from the Taylor formula and (|37|l in Appendb WI Cl 

d 



dL. 



■A,(L 3 ) 



f 1 d 2 

/ dtt—^g(p,m,tL 3 + R p . s ] 

JO ÜL 3 



<c 3 p 5/2 ln(l+p/2) 



Hence, if condition i|31|l is satisfied, then |^|-L2(£3,02) — 1| < 1 and so the lst equation in l|3U|) can be solved to 
express L2 as a differentiable function of L3 and 82- 



C. About Gauss sums. 



Let m, p be relatively prime integers, z an arbitrary complex number, and 

p 

P(p, m,z) = J2 C(p, m, n) z n ; C(p, m, n) = e*™^»- 1 )/? ; Ps = e ^m[2 S+x (p)]/ P _ (32) 

n=l 

where x(p) = 1 when p is even, %(p) = when p is odd, and s = 0, 1, . . . , p — 1. Replacing z — é L in the polynomial 
P(z) one obtains the Gauss sums Q(p, m, L) in Q. The phases of the p s are just the resonant vàlues enumerated 
in a different way In this Appendix we derive the following elementary properties: 

P(p,m,p B+1 ) = p^Pfax^p.) , (33) 

P'(p,m,p ) = i(p + l)P(p,m,p ) foroddp, (34) 

P'(p,m,p ) = — p[P(p,m,p ) + l] forevenp, (35) 
2/00 

\P(p,m, Ps )\ = y/p. (36) 

In addition we derive the following estimates, vàlid for arbitrary z with \z\ = 1: 

|P'(p,m,z)| < Cl p 3 / 2 ln(l + p/2) , \P"{p,m,z)\ < c 2 p 5 / 2 ln(l + p/2) , (37) 

for suitable numerical constants c\, C2, where primes denote derivatives with respect to z. No attempt is made here to 
optimize the bounds l|37|l . and the logarithmic corrections are likely to be artifacts of our proof. 133l) ... l|37|l translate 
in obvious ways into results for the Gauss sums Ç(p,m, L), and their derivatives with respect to L, which were used 
at various places in the main text. Throughout the following we denote w = e l27Tm ^ p , so p s = w s pq. The integers p,m 
being fixed once and for all, we omit specifying them in the arguments of P(.) and C(.). 
Proof of if35|) ,P^P, and if%5|) : from the dehnitions in it is clear that 

C(p-n + l) = (-l) p+1 C{n) , C(n- 1) = w 1 ^ C{n) . (38) 

The first of these identities immediately yields 

P(z) = (-zf^Piz- 1 ) , (39) 

and the second identity yields 

P(zw) = z^Piz) - 1 + z p (-l) p+1 , (40) 

as may be seen from 

p-i p 
P{z) = C(n)z n - 1 + z p (-l) p+1 = ^C{n- l)^ 1 - 1 + z p (-l) p+1 . (41) 

Tl — 71—1 
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Ean. (|40|l in particular yields (|33() . Differentiating (|39|) in z — 1 we obtain 

P'(ï) = Pi 1 ) for odd P > ( 42 ) 
which immediately yields 1|34|1 . because po = 1 when p is an odd number. From 139(1 and 140JI . 

P{zw) = z^-iy+^Piz- 1 ) + 1] - 1 

whence, replacing z — ziw^ 1 ^ 2 : 

P{z 1 w 1 ' 2 ) = z^-lf+^P^w 1 ' 2 ) + 1] - 1 . 
Differentiating in z\ — 1 we obtain: 

PW/2) = — P [Píw 1 / 2 ) + 1] for even p . (43) 
which yields l|35(l because w 1 / 2 = po whenever p is even. 

In order to prové H37|) we need an estimate concerning arbitrary complex polynomials of the form Q{z) = q s z s . Let 
a s = r y s ao, where «o is an arbitrary complex number with |oo | = 1, and 7 = e 27rí /P; anc [ denote Qo = max s |Q(a s )|. 
Then, for any z with \z\ = 1, 

Qo^Q'iz)] < 1 + ip{C + ln(iV + 1)} , (44) 

where C — 3.39968... and N is the integer part of p/2. This may be proven as follows. If r is an integer so that 
-p < r < p then a s = <K r )P: so 

1 P 

g s = -^g(a r )a; s . (45) 

whence the "interpolation formula" follows: 

Q{*) = l E Q{a s )F{za~ 1 ) , F(z) = £ z n = z(z? l)(z l)" 1 . (46) 

" s=l n=l 

If |z| = 1 and z ^ a s for any s, then we denote a the one of the a s that precedes z on the unit circle oriented 
counterclockwisc. Taking derivatives in (|46|l . we obtain: 

n P n N+l 

\Q'(z)\ < ^E^W)!^ E l^«- 1 7- r )l, (47) 
" s=l P r=—N 



where N is the integer part of p/2. Noting that |F'(z)| < p(p + l)/2, 

P 



fiV+l -1 i 

Qo^lQ'^l < p + l + p- 1 E+ E 

L r=2 r=-NJ 
N 

< 



p+i+p- 1 ^! 



|z — Qí7 r | |z — Oí"f r \ 2 

1 



— ' [sin(7rr/p) sin (irr/p) 



r=l 



which directly leads to with C = 2 + i? + 7r 2 /12 where £ = 0.577... is Euler's constant. □ 

Proof of \3f))) : with Q(z) = -P(z), and ao = po, l|33(l shows that |Q(a s )| is independent of s. On the other hand, 
PYX \ls\ 2 = J2i \Q( a s)\ 2 follows from Then JSEl in turn follows, because \q s \ = 1 when Q(z) = P(z). 

Proof of the lst bound in OTp : choosing Q(z) = P(z), yields Q = p 1 ^ 2 and then l(3"7j) follows from (|4"4"ll . 
Proof of the 2nd bound in Tajty : taking the 2nd derivative with respect to z in 146|) . the 2nd derivative of the function 
F(z) appears on the rhs of l|47|) in place of the lst one. Proceeding in a similar way as in the proof of (|44|l . an estimate 
for |Q"(2;)| is obtained, which, using l|36|l leads to l|37fl . □ 
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D. Estimating the Size of an Island. 

We denote 8*, 8* the stable and unstablc cquilibrium positions of 1|12|) in [0,27r]. The separatrix motion occurs at 
energy eV(8*) between the point 8* and the return point 8*, which is the solution of V(8*) = V(8*) with 8* S [0, 2n), 
8* ^ 8*. This orbit attains its maximal momentum at 8 — 8*, so its maximal excursions in momentum and position 
are respectively given by: 

SL 3 = 2y/2p-i\e\[VW)-V(e* s )] , 68 3 = \8* - 8* r \; . (49) 

Introducing a parameter A as in the main text, one may write (|49|l as 

ÍL 3 = 4p- 1/4 |*| 1/a V^(A) ( 68 3 =u(X), 
where the function h(X) is defined as in 

h(X) = A(arcsin(A) - n/2) + \J\ - A 2 (50) 
the value of arcsin being taken in [0, n/2], and u(X) is the continuous function that is implicitly defined by 

Xu = Asin(w) + v 7 ! - A 2 (l - cos(u)) . (51) 

The area of an island is then estimated by A ~ c5L 3 S8 3 with c a slowly varying factor of order unity. This yields lltj|) 
upon defining f(X) = 4u(A) J h(X) . The asymptotics (|17|) in turn follow from the above definitions of u(A) and h(X). 
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This denotation is meant to emphasize that the limit affording this description in terms of trajectòries of a classical 
dynamical system is not the classical limit proper, h — * 0. 
[23] The convention we use in this paper for the signs of Q, k, J is different from those that were adopted in papers appeared 
so far. The positive momentum direction is here the direction of gravity when e > 0, and is the opposite direction when 
e < 0. Though artificial on physical grounds, this choice allows to use the single map Q both for e > and for e < 0. As 
a consequence, in this paper k shares the sign of e, and Í2 and m are always positive. The map that is obtained from Q 
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by changing the sign of k is conjugate to Q under 9 — > 6 + n, so the periòdic orbits of either map one-to-one correspond 
to the periòdic orbits of the other. Thus the tongues are invariant under k — > —k, and this is why \k\ and not k is shown 
on the vertical axis. 

[24] provided that acceleration be always assumed negative in the direction of gravity. 

[25] The jumping index has the physical meaning of a momentum increment. Therefore, if the momentum is assumed positive 
in a fixed direction, independent of the sign of e, then the index is consistently written sgn(e)m when the chosen direction 
is that of gravity and sgn(— e)m in the opposite case. 

[26] The numerical calculations of quantum accelerator modes mentioned in this section (including those in Fig[SJ consist of 
simulations of the exact quantum dynamics of the àtoms, and not of the e-classical dynamics. 

[27] More conventional denotations refer to "branches" in the "Farey tree" . 

[28] A number may be repeated several times in the sequence of Farey approximants constructed that way. 
[29] As Q = GT 2 g/(2ir) is constant along such an EP, the gravity acceleration g has to vary with e. Experimental techniques 
of creating a variable artificial gravity have been devised|2(i|. 



